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Abstract 

The first-order general relativistic theory of a generic dissipative (heat- 
conducting, viscous, particle-creating) fluid is rediscussed from a unified co- 
variant frame-independent point of view. By generalizing some previous works 
in the literature, we derive a formula for the temperature variation rate, which 
is valid both in Eckart's (particle) and in the Landau-Lifshitz (energy) frames. 
Particular attention is paid to the case of gravitational particle creation and 
its possible cross-effect with the bulk viscosity mechanism. 
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I. INTRODUCTION 



The formulation of the cosmological problem, for any particular model, is based on 
a set of phenomenological quantities such as energy density, pressure and temperature, 
which in principle, need to be defined with respect to local observers. In this regard, as 
pointed out by Geheniau and collaborators M, the interdisciplinary field of cosmology is 
the natural inheritor of the concepts and methods ordinarily employed in hydrodynamics 
or, more generally, in non-equilibrium thermodynamics @||. 

Besides being essential to complete the aforementioned phenomenological description, 
the thermodynamic analysis of cosmological models may also have interesting observational 
consequences. In principle, the dimensionless entropy content of the (observable) universe, 
a huge number of the order of 10 87 , may have been generated by dissipative mechanisms 
like heat conduction, viscosity, "chemical reactions", diffusion, photon creation, etc., taking 
place in an inhomogeneous and anisotropic stage of the early universe [|]||. In this case, after 
some few minutes, the gradients of the main physical quantities will be constrained from 
primordial nucleosynthesis studies. Later on, at the time of recombination, these irreversible 
mechanisms may also contribute to the temperature anisotropy of the cosmic background 
radiation. Indeed, although severely restricted from COBE measurements, the presence 
of the classical dissipative mechanism are naturally expected in a non-homogeneous quasi 
FRW universe, and, in principle, their specific signature must be somewhat different of the 
standard Sunyaev-Zeldovich and Sachs- Wolf effects ||. 

In order to unify in a single and coherent scheme all irreversible phenomena occurring 
in a simple fluid or in mixtures, standard non-equilibrium thermodynamics works with two 
basic ideas @||. The first one is the local equilibrium hypothesis, whose mathematical 
expression is given by the equilibrium Gibbs law in its local form. It implies that, out 
of equilibrium, the basic state functions such as the entropy, depend locally on the same 
set of thermodynamic variables as in equilibrium. In particular, the usual thermodynamic 
temperature and pressure concepts are maintained in the non-equilibrium regime. The 
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second idea is that, in the presence of dissipative processes, there is a local entropy source 
strength r (entropy variation per unit volume and unit time), which, by the second law of 
thermodynamics, is always non-negative. Mathematically, it takes the form of a balance 
equation with r as a source term. 

By combining such assumptions with the fluid equations of motion, one finds an ex- 
pression for the entropy source strength, as well as for the constitutive (phenomenological) 
relations themselves. The usual equilibrium theory is readily recovered by taking the limit 
of vanishing entropy production rate. 

This approach has systematically been applied at the classical level |||| as well as in 
the special-relativistic domain ||- ||. Its extension to the general-relativistic framework 
is straightforward, provided the gravitational field varies slowly over the mean free path or 
the mean free time of the fluid particles |4j]. Formally, since, under these conditions, the 
equivalence principle should hold, the generally covariant equations may be established by 
the usual minimal coupling recipe, that is, the replacement of usual derivatives by covariant 
derivatives and the replacement of the Minkowski metric r] a p by its Riemannian counterpart 
g a/ 3. However, as is widely known, in the relativistic case, there is an ambiguity related to the 
possible choices of the macroscopic hydrodynamic four- velocity. In Eckart's formulation [|7|], 
the four-velocity is directly related to the particle flux, while, in Landau-Lifshitz's approach 
||, it is directly related to the energy flux. In principle, a general treatment should be able 
to deal with any of these "gauge" choices. 

In this paper, we consider a unified covariant description of a heat-conducting viscous 
simple fluid with particle creation. The theories of Eckart and Landau-Lifshitz will be seen 
to be special cases of this general formulation. We also derive a new temperature evolution 
law which holds for any particular choice of the macroscopic hydrodynamic four-velocity. 
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II. GENERAL RELATIVISTIC THEORY OF FLUIDS 



The thermodynamic state of a relativistic simple fluid is characterized by an energy- 
momentum tensor T Q/3 , a particle current N a and an entropy current S a . The fundamental 
equations of motion are expressed by the conservation law (semi-colon denotes covariant 



derivative) of energy-momentum and the equation of balance for the particle number |L0 

1*00 = , (1) 

N a ; a = V , (2) 

where \I/ is a particle source (\I/ > 0) or sink < 0) term. The second law of thermody- 
namics requires that the entropy source strength be non-negative 

S a ;a = T>0 , (3) 

where r = describes a non-dissipative state, and r > denotes a dissipative state. A per- 
fect fluid always evolves through non-dissipative (equilibrium) states, whereas an imperfect 
fluid typically evolves through dissipative (non-equilibrium) states (see, however, |TT1| ). 



A. Adiabatic Limit 

By choosing an arbitrary hydrodynamic frame of reference, whose four-velocity obeys 
u a u a = 1, the primary variables T a/3 , N a and S a take the following forms |5|,[|] 

T a ? = pu a uP -ph a ?, (4) 
N a = nu a , (5) 

S a = nau a , (6) 

where the tensor := u a, uP — g af3 is the usual projector onto the local rest space of u a . 
The variables p, p, n and a stand respectively for the energy density, thermostatic pressure, 
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particle number density and specific entropy (per particle), and are related by the so-called 
Gibbs law §g 

nTda = dp — - - dn . (7) 

n 

As a particular result of the general theory, imply that the entropy source strength 

defined by (|3]) vanishes identically, 

S a ;a = r = , (8) 
as should be expected for thermal equilibrium states. 

B. Non-Equilibrium States 

In principle, the inclusion of dissipative processes such as heat conduction and viscosity, 
requires additional terms in the primary variables describing a perfect fluid. However, unlike 
in the adiabatic limit case, the presence of a heat transfer poses a problem regarding the 
definition of the hydrodynamic four-velocity u a . It is necessary to specify whether u a is the 
four- velocity of the energy transport or particle transport. In Eckart's formulation, u a is 
identified with the four- velocity of particle transport (particle frame) j7|. In the approach 
of Landau-Lifshitz, u a is defined as the four- velocity of energy transport (energy frame) || . 
Formally, the particle frame is the unique unit time-like vector parallel to N a , whereas the 
energy frame is the unique unit time-like eigenvector of T af3 . Both theories assume that, for 
weak space-time gradients, the basic quantities contain no terms higher than first order in 
deviations from equilibrium. 

In the presence of irreversible processes, we must add small terms AT a/3 and AN a in (|J) 
and (|^), which are restricted by the second law of thermodynamics (|j) 

T a/3 = pu a vP - ph a(3 + AT al3 , (9) 
N a = nu a + AN a , (10) 
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Usually, at this point, one specifies whether the Eckart or Landau-Lifshitz approach will 
be adopted. However, this is not necessary in the covariant frame-independent formulation 
presented here since these theories will be recovered as particular cases. 

As before the fluid motion equations are contained in (|1|) and (|2|), which express the 
energy conservation law and the balance equation for the particle number, respectively. 
Therefore, differentiating the expression and projecting it on the direction of the four- 
velocity u a , one finds 

u a T^. p = p+(p + p)9 + u a AT^. p = , (11) 

and from (|2]) 

N a . a = h + n9 + AN a . a = m , (12) 

where an overdot means the derivative along the world lines of the fluid volume element, 
e.g., p = uPp-fi, and 9 — u^.p is the expansion rate of the fluid. 

Now, taking the covariant derivative of (0) along the world lines of the fluid volume 
element, and making use of equations ( |TTD and (0), we obtain 

T{nau%!3 = -« a Ar% - pA> + /iAiV% , (13) 

where /i is chemical potencial defined by Euler's relation 

li =E±P-oT . (14) 
n 



Defining the entropy flux 



S^nau^-^AN^ + ^AT^ , (15) 
we see from (|T3|) that the entropy source strength assumes the following form 

S%= (^_Z|^)A^-(f _^)AiV»-^ , (16) 
which is a function dependent only on the dissipative fluxes. 



At the level of the primary fluxes, the effect of the dissipative processes is to add up new 
fields for the energy-momentum tensor and the particle flux vector [fT||5|, |T0| , |13|| , which must 
be constrained by the second law. We have 

AT a ? = -(Tl + p c )h a P + q a U P + qPu a + Tl a P , (17) 

and 

AN? = z/ , (18) 

where the five additional fields H,p c , q a , and IT a/3 are, respectively, the bulk viscous pres- 
sure, the creation pressure (due to the gravitational matter creation |10|), the heat flow, 
the particle drift, and the shear viscosity stress. These fields describe the deviations from 
equilibrium within the fluid and satisfy the following constraints 

u a q a = u a v a = u a U a(] = g a ^ = IL™ = , (19) 

where the square brackets denote anti-symmetrization and round brackets, below, sym- 
metrization. 

As one may check, inserting equations ( |T7|) and ([18]) into (|T6|), and using (|19D , it is readily 



seen that 



_ m p c e (T.p up\n (hp ij.t. p \ p u a ,^ 

o - a --— ~Y ~r~ VT2 Y) VT ~T r ) T { ' 



On the other hand, the covariant derivative of u a may be decomposed as [14 



u a -p = a a p + ^9h a p + ii a up + u af3 , (21) 



where 



^ = \Kh^u^, u) -hh, v ) , (22) 



and 



u a p = -h^h" p u M , (23) 
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are respectively the symmetric traceless shear tensor and the vorticity tensor. Inserting (]2T 
into equation (|20|), we may rewrite the entropy source as 



not ^ \^;p ™pj r\p h"^;p \ H , — » ap 

S ^ = T S — ^— - —) V P + - , (24) 

where t s describes the entropy source strength due to the scalar irreversible processes 

r> = _™_*»_£ . (25 ) 

Hence, to the first-order of approximation, ( p4|) will be consistent with the second law of 
thermodynamics if the phenomenological relations among the vector and tensor dissipative 
fluxes and thermodynamic forces are taken to be 

q a = X r ; r = h^(T.p-Tu p ) , (26) 

= ^ (27) 

u a = (X a ; X a = h al3 (^) , (28) 



T 
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where x> V> an d C stand, respectively, for thermal conductivity, shear viscosity and "difusion" 
coefficients. Using (|26|)-(^8|), the entropy source strength ( p4|) becomes 

qa _ ^ _ x^r _ ckx« , w^ afi (9Q , 

and since the heat flow and the particle drift are space-like vectors (0 Q Q < 0, X a X a < 0), 
the second law of thermodynamics will be satisfied if, for any configuration of the fluid, x, 
r\ and £ are positive. In order to treat the scalar irreversible processes, we observe that they 
are fluxes with same tensor rank (scalars), and in the first-order of approximation this give 
rise to cross effects In order to describe this cross effect, we propose the following 



relations among fluxes and thermodynamic forces, 

n=-£n(!)-£ia0 , (30) 



and 

Pc = -£210-62 (|) , (31) 

where II, p c are fluxes and (j\ @ are thermodynamic forces, with £y > (i,j = 1,2) 
phenomenological coefficients. When the particle source strength \1/ is zero, we see from 
(|30l)- (|3~l|) that II = — £ 12 # and p c = —£21$- Recalling that the phenomenological coefficients 
are supposed to obey Onsager's reciprocity relations, it follows that £12 = £21- This means 
that in this particular case (no particle creation) the creation pressure and bulk viscosity 
reduce to the same process. In the general case, the entropy source (|29|) can be written as 



& ia — rp 2 + rp rp 2 + rp + rp \ 6Z >) 

We see that the second law of thermodynamics will be satisfied for any configurations of the 
fluid if XiViC£i2 are positive definite. We also see that the remaining coefficients £n and £22 
must satisfy the following inequalities: £n + £22 > ^ if ^ > 0, and £n + £22 < if ^ < 

(cf. H). 

Concerning the results presented in the reference we have the following remarks. 
First, we observe that the flux associated with the particle creation rate, namely, the creation 
pressure was not considered. Second, the authors have used II as a force and not as, we think, 
a thermodynamic flux. Another question is related to the hypothesis implicitly adopted in 
their paper, that the creation pressure and the bulk viscous pressure are the same (II = 
p c ). In particular, from the phenomenological laws obtained there, we see that for photons 
(/i = 0), the creation pressure is given by p c = however for the static fluid 9 = 0, 

and using the energy conservation, we find that p = 0. These results seem to be physically 
inconsistent, because we have energy density p constant and simultaneously photon creation, 
which should be responsible for a variation of p. Note also that, in the adiabatic limit, the 
entropy source strength vanishes, as expected (see §). 

Summing up, the basic set of equations governing the first-order theory of dissipative 
simple fluids in a general formulation are: the energy conservation law ([]]) and equation of 
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balance (g) with T af3 and N a given by (|9|) and (|10|); the Gibbs law (0), the constitutive 
equations (|26|)-(|28|) and the constitutive equation for a scalar process (|30"|)-(|31|) with the 
fluxes constrained by (|1^) and linked with the entropy source strength by (0). Notice that 
such a system is underdetermined since, in principle, there are fifteen independent equations 
and seventeen unknowns; more specifically, as in the classical theory, the time behavior of 
the relativistic one-component fluid can only be determined, for specified initial boundary 
conditions, by adding two equations of state, say 

p = p{n,T) , (33) 

and 

p = p(n,T) . (34) 

We notice that, although sucessful in revealing the physics underlying a large class of phe- 
nomena, the first-order theories present some experimental and theoretical drawbacks. In its 
classical version, the linear constitutive equations (|2B|)-(PBD and (^Dj)-(^) are not adequate 



at high frequencies or short wave lengths as manifested in experiments on ultrasound propa- 
gation in rarefied gases and on neutron scattering in liquids |T(| . Besides they also allow the 
propagation of perturbations with arbitrarily high speeds, which, although perhaps merely 
unsatisfactory on classical grounds, is completely unacceptable from a relativistic point of 
view; furthermore, they do not have a well-posed Cauchy problem and their equilibrium 
states are not stable. Several authors have formulated relativistic second-order theories 
which circumvent these defficiencies P,[l7|-|T9|1 . In a forthcoming paper [20|, we intend to 



extend our considerations to this class of theories. 



III. TEMPERATURE EVOLUTION LAW 

Let us now derive the general equation for the temperature law in a relativistic fluid 
taking into account all the irreversible processes. In what follows we consider that the 
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energy density and the equilibium pressure are functions of the thermodynamic variables n 
and T. Differentiating p = p(n, T) along the world lines of the fluid volume element one 
finds 



(35) 



and combining with the energy conservation law ( JTT| ) and the balance equation (|T2P, we 
obtain 



v QJ>> n 



(|£)t(* - AiV%) - u a AT°' , p 



nap 



(36) 



Since da must be an exact differential, (|7|) leads to the thermodynamic relation 

'p + p" 



d_ fj_ 
df [nT 



dp 
dn 



T 
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or equivalently, 



d_ fj_ (dp_ 
dn~ [nT [df 



dp s 



n/ AT 



(37) 



(38) 



and combining 
temperature 

t _ 
f ~ 



with (|36|), we obtain the general equation governing the variation of 



^ (* - AiV%) + u a AT^, p 



(39) 



where all dissipatives fluxes are described by AA^ a and AT af3 . In the adiabatic limit, the 
above equation reduces to the temperature variation rate of a perfect fluid 



T 
f 



' dp\ 



(40) 



Consider now in the Landau-Lifshitz or energy frame. In this case, the comoving 
observers do not see the irreversible contribution of the energy flux, that is, q a = in the 
energy- momentum tensor. As we have seen, there is an additional contribution v a in the 
particle flux AN a . We obtain 

dp, 



r r)T ' n 



dn 



)t« - *) - (n + p c )6 LL + a af3 Tl a P 



(41) 
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where 

— n — v a 

e LL = ;Q (42) 

n 

is the expansion of the fluid in the Landau-Lifshitz frame. 

Now let us consider the Eckart or particle frame formulation. In this case, the comoving 
observers do not see the irreversible contribution of particle drift, that is AN a = v a = in 
the particle flux vector. We have 



)vPE + 



T ~ y dp ,n * T($) 
where 



-(p c + U)9 E + ii a q a + a aP U^ - q?. fi - (^) T * 



(43) 



Be = (44) 
n 

is the expansion of the fluid in the Eckart frame. 

It should be noticed that in the adiabatic limit, when all dissipative fluxes are absent, 
equations ( f£Tp and ( pE3"D reproduce the same temperature variation rate of a perfect fluid. The 
same happens when the dissipative fluxes reduce to a creation pressure plus the bulk viscosity 
mechanism. For this homogeneous and isotropic dissipative simple fluid, the temperature law 



has previously been derived by Lima and Germano (see also Zimdahl |22| to the case of a 
two-fluid mixture). We also emphasize that, for each formulation, the complete temperature 
evolution equation must play an important role in problems of astrophysics and cosmological 
interest envolving the classical non-equilibrium mechanisms. In the cosmological domain, for 
instance, a natural application of equations flilD and fl43|), is related to the contributions of 



the irreversible processes to the anisotropy temperature of the cosmic background radiation 
during the decoupling time. 
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